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, Abstract 



superprocess limit for an interacting birth-death particle system modelling a population 
^ ' with trait and physical age-structures is established. Traits of newborn offspring are inherited 

from the parents except when mutations occur, while ages are set to zero. Because of 
I interactions between individuals, standard approaches based on the Laplace transform do 

not hold. We use a martingale problem approach and a separation of the slow (trait) and 
fast (age) scales. While the trait marginals converge in a pathwise sense to a superprocess, 
the age distributions, on another time scale, average to equilibria that depend on traits. The 
convergence of the whole process depending on trait and age, only holds for finite-dimensional 
time-marginals. We apply our results to the study of examples illustrating different cases of 
' trade-off between competition and senescence. 
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\^ ■ 1 Introduction 



We consider an asexual population in which the survival probability and reproduction rate of 
each individual are characterized by a quantitative trait, such as for example the body size, or 
the rate of food intake. As emphasized by Charlesworth [7], most of these abilities also depend 
^ ' on age. In this paper, we are interested in studying the joint effects of age and trait structures 

in the interplay between ecology and evolution. Evolution, acting on the distribution of traits in 
the population, is the consequence of three basic mechanisms: heredity, which transmits traits 
1^ ■ to new offspring, mutation, which drives the variation in the trait values, and selection between 

j_j ■ these different trait values, which is due to ecological interactions. Some questions on evolution 



are strongly related to the age structure. For example, we would like to understand how the age 
influences the persistence of the population or the trait's evolution, or which age structure will 
appear in long time scales for a given trait. 

Our model relies on an individual-based birth and death process with age and trait introduced 
in Meleard and Tran [26] (see also Ferriere and Tran [14j). It generalizes the trait-structured case 
developed in Champagnat et al. [5l[6] and the age-structured case in Jagers and Klebaner [TSlfTQ]. 
Tran [34J. Here, each individual is characterized by a quantitative trait and by its physical age, 
that is the time since its birth. We describe the dynamics of large populations composed of small 
individuals with short lives. Life-lengths and durations between reproductions are assumed to 
be proportional to the individuals' weights, and these weights are inversely proportional to the 
population size. At each birth, the age is reset to 0, inducing a large asymmetry between 
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the mother and her daughters. When a mutation occurs, the new mutant trait is close to its 
ancestor's one, yielding a slow variation of the trait. Hence, these two mechanisms lead to a 
difference of time scales between age and trait. Moreover, we take resource constraints into 
account by including competition between individuals. 

Our aim is to study an approximation of this model when the population size tends to 
infinity. Our main result shows that two qualitatively different asymptotic behaviors arise from 
the separation of the fast age and slow trait time scales. While the trait marginals converge in 
a pathwise sense to a superprocess, the age distributions stabilize into deterministic equilibria 
that depend on the traits. To our knowledge, nothing has been done before in the setting we 
are interested in, with slow-fast variables and nonlinearity. The techniques we use are based 
on martingale properties and generalize to this infinite dimensional setting the treatment of the 
slow-fast scales for diffusion processes, developed by Kurtz [23], and by Ball et al. [2]. 

Our results generalize Athreya et al. [1], Bose and Kaj [3], where averaging phenomena 
are proved in the case where birth and death rates do not depend on age. In our case with 
dependence, the lifelength of an individual cannot be governed by an age distribution function 
independent of trait, or by a positive continuous additive functional, as in Bose and Kaj [3], 
Dawson et al. [9], Dynkin [lOj, Kaj and Sagitov jj^, Fleischmann et al. [16j, Wang [35]. More- 
over the Laplace characterization that these authors use extensively does not hold anymore when 
interactions between particles are allowed. In Evans and Steinsaltz |13j . the damage segregation 
at cell fissioning is considered ge. Nevertheless in their model there is no interaction 

between cells, and at each birth, the daughters' ages (as damages) are not reset to zero, but 
distributed asymmetrically following a distribution centered on the mother's age. 

In Section [21 the population dynamics is described by an individual-based point measure- 
valued birth and death process with an age and trait structures. In Section O we establish our 
main limit theorem (Theorem 13. ip where the averaging phenomenon is obtained. In Section 21 
we consider as an illustration two models where the population is structured by size and physical 
or biological age. We present simulations and comment the different behaviors. 

Notation: For a given metric space E, we denote by B([0,T],£') the space of right continuous 
and left limited (cad-lag) functions from [0, T] to E. This space is embedded with the Skorohod 
topology {e.g. [Ml ED]). 

If is a subset of M*^, we denote by Jv[p[X) the set of finite measures on X, which will be 
usually embedded with the topology of weak convergence. Nevertheless, if X is unbounded, we 
will also consider the topology of vague convergence. If we need to differentiate both topological 
spaces, we will denote by (A^i?(Af), w), respectively [M.p{X).,v)., the space of measures endowed 
by the weak (resp. vague) topology. For a measurable real bounded function /, and a measure 
H G A^ir(Af), we will denote 

(/^,/)= / f{x)fi{dx). 
Jx 

For £ G N, we denote by C^{X, M) the space of real bounded functions / of class with bounded 
derivatives. In the sequel, the space q'^{X X M+,M) (resp. CciX,R), Cl{R+,R)) denotes the 
space of continuous bounded real functions ^p{x, a) on X x of class with respect to a with 
bounded derivatives (resp. of continuous real functions on X with compact support, of real 
functions on M_|_ with compact support in [0, -|-oo)). 
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2 Microscopic age and trait structured particle system 



We consider a discrete population in continuous time wliere the individuals reproduce, age and 
die with rates depending on a hereditary trait and on their age. An individual is characterized 
by a quantitative trait x & X where X is a closed subset of M*^ and by its physical age a G M_)_, 
i.e. the time since its birth. The individuals reproduce asexually during their lives, and the trait 
from the parent is transmitted to its offspring except when a mutation occurs. Resources are 
shared by the individuals, implying an interaction described by a kernel comparing the competi- 
tors' traits and ages. Senescence, which quantifies the decrease of fertility or survival probability 
with age, is also taken into account. These two phenomena create selection pressure. 

We are interested in approximating the dynamics of a large population whose size is parametrized 
by some integer n. This parameter can be seen as the order of the carrying capacity, when the 
total amount of resources is assumed to be fixed. If the parameter n is large, there will be many 
individuals with little per capita resource and we renormalize the individual biomass by the 
weight 1/n. 

We consider here allometric demographies where the lifetime and gestation length of each in- 
dividual are proportional to its biomass. Thus the birth and death rates are of order n, while 
preserving the demographic balance. As a consequence the right scale to observe a nontrivial 
limit in the age structure, as n increases, is of order 1/n. 

The population at time t is represented by a point measure as follows: 



where = {nX^, 1) is the number of individuals alive at time t, and Xi{t) and Ai{t) denote 
respectively the trait and age of individual i at time t (individuals are ranked in lexicographical 
order for instance). 

The dynamics of X"' is given as follows: 

• The birth of an individual with trait x E X and age a G M+ is given by nr{x, a) + b{x, a). 
The new offspring is of age at birth. Moreover, it inherits of the trait x of its ancestor with 
probability 1 —p{x, a) G [0, 1] and is a mutant with probability p{x, a) G [0, 1]. The mutant 
trait is then x + h, where the variation h is randomly chosen following the distribution 



• Individuals age with velocity n, so that the physical age at time t of an individual born at 
time c is a = n{t — c). 

• The intrinsic death rate of an individual with trait x e X and age a G M+ is given by 
nr{x,a) + d{x,a). The competition between individuals {x,a) and is described 
by the value C/((x, a), (y, a)) of a kernel U. In a population described by the measure 
X & X the total interaction on an individual {x,a) is thus: 




(2.1) 



1=1 



7r"(x,d/i). 




(2.2) 



and its total death rate is nr{x, a) + d{x, a) + XU{x, a). 



Assumption 2.1. 1. The birth and death rates b and d are continuous on X x and 
bounded respectively by b and d. 
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2. The function r is continuous on X x M+. There exist a positive constant r and a non- 
negative real function r such that V(x,a) G A" x R_|_, r(a) < |r(x,a)| < f with 



+00 



r{a)da = +00. 



(2.3) 



3. The competition kernel U is continuous on {X x and is bounded by U. □ 

Assumption 12. II - (2) implies that any individual i from the population X" has a finite lifetime 
that is stochastically upper-bounded by a random variable Df{t) with survival function 

pi 

S'^ii) = F{D^{t) > £) = exp f - / nr{nu)du) , (2.4) 



where we recall that the aging velocity is equal to n. 

If the competition kernel U is positive on {X x M+)^, it can model a competition of the logis- 
tic type: the more important the size of the population is and the higher the death rate by 
competition is. For examples of such kernels we refer to [26j. 

Example 2.2. Let us illustrate the condition (j2.3|) . 

1. If the function r is lower bounded by a positive constant r, then (j2.3|) is satisfied and so is 
()2.4p . with exponential random variable D^{t) of parameter nr. 



Another example is when the trait x is linked to the rate of metabolism, which measures 
the energy expended by individuals. Ageing may result from toxic by-products of the 
metabolism and we can define a biological age, xa. If x G [xi,X2] with xi,X2 > and if 
we define r(x,a) = xa, then Condition (|2.3p is satisfied with r(a) = xia. The example of 
Section 14.21 deals with this case. 



3. If we consider r(x,a) = 7/(1 -|- a) with 7 G (0,1), then ()2.3p is also satisfied and the 
probability of observing an age higher than a is equivalent to when a tends to infinity. 
Such cases with distributions in the domain of attraction of a stable law, but without 
interaction, have been considered for instance in |16j . 

Assumption 2.3. For any x £ X, the mutation kernel 7r^{x, dh) has its support in ^ — {x} = 
{/i G M°' I X -|- /i G X}. We consider two cases: 

1. The trait space X is a compact subset of and there exists a generator ^ of a Feller 
semi-group on Cb{X,M.) with domain 'D{A) dense in Ci,{X,U.) such that: 



V/ G V{A), lim sup 



n I {fix + h)- /(x))7r"(x, dh) - Af{x) 

X-{x} 



0. (2.5) 



2. The trait space X is a. closed subset of and we assume in addition that there exists 



?i > 4 > 2 with qi(A',R) C V{A) and V/ G (Af,M), Vx G X. 

\Af{x)\<C 



(2.6) 



and 



sup 



|fc|<£o 
k={ki,...,ka) 



n I {fix + h)- /(x))7r"(x, dh) - Af{x) 
X-{x} 



(2.7) 



|fc|<£i 
fc=(fci,...,fcd) 



where D^f{x) = dl\ . . . d^jf{x), En is a sequence tending to as n tends to infinity and 
C is a constant. 
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Remark 2.4. Both Assumptions ()2.5p and (j2.7p describe small mutation steps. The stronger 
hypothesis ()2.7p is required when X is not compact, to obtain the tightness in the proof of 
Theorem 13.11 

Example 2.5. Let us give some examples of mutation kernels satisfying ()2.5p or ()2.6p and ()2.7p . 

1. In the case where = [xi, X2], the mutation kernel 7r"(x, dh) is a Gaussian distribution with 
mean and variance cr^/ra, conditioned to [xi — x, X2 — x]. In this case, elementary computation 
shows that for / G C^{[xi,X2],R) such that /'(xi) = f'{x2) = 0, = ^f"{x), which 
satisfies (|2.5p . 

2. In the case where X = W^, a possible choice of mutation kernel 7r'^{x, dh) is a Gaussian 
distribution with mean and covariance matrix S(x)/n, with = < i,j < d). 
The generator A is given for / € Cj2(k'^^m) by = ^ I]f,j=i ^ij{x)dfjf{x). If the function 
S is bounded, then Assumption (j2.6p is fulfilled. If moreover, the third moments of 7r'^{x, dh) 
are bounded (in x), then ()2.7p is satisfied. □ 

Let us now describe the generator L" of the AiF{X x ]R_|_)-valued Markov process X", which 
sums the aging phenomenon and the ecological dynamics of the population. As developed in 
Dawson [8] Theorem 3.2.6, the set of cylindrical functions defined for each fi £ Aip{X x by 
F^{n) = F{{fi,ip)), with F £ Cl{R,R) and if £ C°'^{X x M+,M), generates the set of bounded 
measurable functions on Mf{X x M+). For such function, for fi G AiF{X x M+), 



Under the condition sup„gj^, E {{Xq, 1)) < +00, it has been proved in Meleard-Tran [26] (see also 
the case without age in Champagnat-Ferriere-Meleard jB] and the case without trait in Tran [33] ) , 
that there exists for any re, a cad-lag Markov process with generator L", which can be obtained 
as solution of a stochastic differential equation driven by a Point Poisson measure. Trajectorial 
uniqueness also holds for this equation. The construction provides an exact individual-based 
simulation algorithm (see [Hj). 

A slight adaptation of the proofs in [6] allows us to get the 
Proposition 2.6. (i) Under Assumptions \2.1[ and if 



n{fi,daipi.))F'{fi) 




where 



K"(a;, a, dh) = p{x, a) 7r"(x, dh) + (1 - p{x, a)) 6o{dh). 



(2.9) 



sup E((XoM)3 



) < +00 



(2.10) 



neN* 



then for all T > 



sup sup e((X",1)^) < +00 and sup e( sup {X^,lf) < +00. 



(2.11) 
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(n) Moreover, for n G N* and a test function tp G C^'^{X x M_|_,M), the process M"'^ defined by 



Mf''^ = (Xf,v9)-(Xo",(^)-n / {X^,da^{x,a))ds 







/ ([nr{x,a) + b{x,a)) / ip{x + h,0)K^\x,a,dh) 

JxxR+ ^ Jr'^ 



- {nr{x,a) + d{x,a)+Xl'U{x,a))ip{x,a)jX^{dx,da)ds (2.12) 
is a square integrable martingale started at with quadratic variation: 

{M'^'V)t = - [ [ ({nr{x,a) + b{x,a)) [ ip^{x + h,0)K'^{x,a,dh) 

" Jo JxxR+ ^ JR'i 

+ {nr{x,a) +d{x,a) + X^U{x,a))ip^{x,a)^ X^{dx, da) ds. (2.13) 

Notice that in ()2.12p . the mutation rate is hidden in the kernel K": 

(nr(x, a) + b{x, a)) / (p{x + h, 0)K^{x, a, dh) = {nr{x, a) + 6(x, a)) (1 — p(x, a))ip{x, 0) 

+ {nr{x,a) + b{x,a))p{x,a) 92(2; + /i, 0)7r"(x, (i/i). (2.14) 

jRd 

3 Superprocess limit 

We now investigate the hmit when n increases to +00. In the hmit, we obtain a continuum of 
individuals in which the individualities are lost. It is in particular difficult to keep track of the 
age-distribution when n tends to infinity. Because of the non-local branching (a mother of age 
a > gives birth to a daughter of age 0) and because the aging velocity tends to infinity, it is 
impossible to obtain directly the uniform tightness on D(M+,7Wi?(Af x M+)) of the sequence of 
laws of the measure-valued processes {X'^{dx, da))neN* , as it can be observed considering (j2.12p . 
Indeed, assuming that the function (p only depends on a, the term of the form 

/ / nr{x,a) {p{0) - p{a)) Xg{dx,da)ds 

Jo JxxR+ 

cannot be tight if r(x, a) is bounded and X" is tight. Therefore, we will be led to firstly show the 
uniform tightness of the trait marginal of the process X" and then to prove that in the limit, an 
averaging phenomenon appears for the age dynamics. Indeed, this "fast" evolving component 
stabilizes in an equilibrium that depends on the dynamics of the "slow" trait component. 

We generalize to measure- valued processes, averaging techniques of Ball et al. [2], Kurtz 
|23j for diffusion processes. A specificity in our case is that the fast-scaling is related to time, 
since age is involved. In addition, notice that the competition between individuals creates a 
large dependence between the age and trait distributions. At our knowledge this dependence 
has never been investigated before in the literature. 

Let us introduce the marginal X^{dx) of X^{dx,da) defined for any bounded and measurable 
function f on X and for any t £ M+ by 



[ f{x)Xndx)= [ fix)X^{dx,da). (3.1) 

Jx JXxR+ 
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Our main result states the convergence of the sequence (X")„gN* to a nonhnear super-process. 
The nonUnearity remains at the slow time scale in the growth rate, which is preserved in this 
asymptotics. Moreover, fast mutations are compensated by small mutation steps. Fast births 
and deaths provide stochastic fluctuations in the limit. 

Theorem 3.1. Assume Hypotheses \2.1\ and \2.'J[ i2. 1 U\) and assume that there exists Xq G 
X M^-) such that lim„_^+oo -'^o ~ "'^o (A^_f('^ x 1^+))?^); limit being in probability 
for the sake of simplicity. 

For any x G X , let us introduce the age probability density 

Jq exp [ — Jq r[x,a)aa)aa 

Then, for each T > 0, the sequence converges in law in ID)([0,T], {Aip{X),w)) to the 

unique superprocess X G C{[0,T],M.p{X)) such that for any function f £ 'C>{A), 



m/ = {Xt,f)-{XoJ)- / ((f^)ix)Af{x) 

+ - {^{x) + xSj{x))\f{x)^X,{dx)ds (3.3) 
is a square integrable martingale with quadratic variation: 

{Mf)t= [ [ 29{x)f^{x)Xs{dx)ds. (3.4) 
JO Jx 

Here, any tp{x) is defined for a bounded function ij){x^ a) by 

'il){x) = / ip{x,a)m{x,a)da, 



and XtU{x) is given by 

XtU{x) = 1(1 I U{{x,a),{y,a))Tn{y,a)da m{x,a)da\ Xt{dy). 



X 



□ 



Theorem 13.11 states that in the limit, an averaging phenomenon happens and the "fast" age 
component finally submits to the dynamics of the "slow" trait component. Since the fast- 
scaling (involving age) is related to the time, the stable age distribution m(x, a)da given in ()3.2p 
is obtained for each trait x as the long time limit in the age-structured population where all 
coefficients except r(x, a) are zero. 

Before proving this slow-fast limit theorem, let us insist on the main difficulty created by the 
competition mechanism. Indeed, the branching property fails and it impedes the use of Laplace- 
transform techniques, as it had almost systematically been done in the past papers studying 
particle pictures with age-structures. Our model generalizes the age-structure population process 
studied in Athreya et al. [T], Bose and Kaj [3], in which birth and death rates are equal to a 
constant A. In that case, the limiting behaviour of the renormalized critical birth and death 
process appears as a particular case of Theorem 13.11 with m(x, a)da = \e~'^°'da. In Dawson 
et al. |9j, Dynkin [10], and Kaj and Sagitov [22j, the age dependence is modelled through an 
additive functional of the motion process. In that way, the age "accumulates" along the lineage. 

7 



In our case, the age is set to zero at each birth, inducing a renewal phenomenon. The life-length 
does not have a fixed probability distribution anymore, unless there is no interaction. In 
the authors consider a particle system with a different scaling, which favors large reproduction 
events. The limit in this case is not a superprocess anymore but behaves as the solution of 
a McKendrick-Forster equation perturbed by random immigration events created by the large 
rare birth events. 

The proof of Theorem 13.11 is the aim of Section [3l We firstly establish the tightness of the se- 
quence {X"-)n£N* (Section EH]). To identify its limiting values X, we need an intermediary step. 
We consider the measures {Xl^{dx,da)dt)n£N* and show that their limiting values are equal to 
Xt{dx)fh{x,a)dadt. This implies that X is solution of the martingale problem given in (j3.3p . 
()3.4p . Uniqueness in this martingale problem allows us to deduce the convergence of 



3.1 Tightness of (X")„gN. 

In this subsection, we shall prove that: 

Proposition 3.2. Assume Hypotheses \2.1\ and \2.Sl and i2.1U\) . If the sequence of laws of 
{XQ)n£N* is uniformly tight in {A4f{X),w), then the sequence (£(X"'))„gpj* is uniformly tight 
in the space of probability measures on'S}{[0,T],{A4F{X),w)). 

Proof. Recall firstly that for a measurable and bounded function f on X, the process 
M^^f = {X^,f)-{X^J)- f I ((nr(x,a) + 6(x,a)) / f{x + h)K^{x,a,dh) 

Jo JXxR+ ^ Jr^ 

- {nr{x,a) +d{x,a) + XgU{x,a)) f{x)^Xg{dx,da) ds 
= (Xr,/)- W,/)- f [ ([b{x,a)-d{x,a)-X^U{x,a))f{x) 

Jo JXxR+ ^ 

+ {nr{x,a) +b{x,a))p{x,a) / {f{x + h)-f{x))TT''{x,dh)\X'^{dx,da)ds (3.5) 

JR<1 ' 

is a square integrable martingale started at with quadratic variation 

(M''-S)t = - ( I ({nr{x,a)+b{x,a)) [ f'^{x + h)K'^{x,dh) 
Jo JXxR+ ^ JRd 

+ (nr(x,a) +d{x,a) + X^U{x,a)) f'^{xj^ X^{dx,da) ds 

1 



n 



2nr{x,a) + b{x,a) + d{x,a) + X2U{x,a)) f'^{x) 

JXxT ' 



+ {nr{x,a) + b{x,a))p{x,a) / {f'^{x + h)-f{x))7r'^{x,dh)]x'^{dx,da)ds. 

(3.6) 

We divide the proof into several steps. 

Step 1 Firstly, WG prove the uiiiforni tightness of (/^(^^))^^ in the space of probability 
measures on B([0,T], (A^i;'(Af),t>)). 

Let us consider a continuous bounded function / E ^(^), and show the uniform tightness of 
the sequence {{X^^ , f))neN* in B([0, T], M). We remark that for every fixed t € [0, T] and n > 0, 



„_,^E(sup^prnri(^r,l)) 
P(KXr,/)| >k)< "^"^ ^ ^^ (3.7) 



which tends to as A; tends to infinity (cf. (|2.1ip ). This proves the tightness of the family of 
time-marginals {{Xf, f))neN*- Denoting by A"'-^ the finite variation process in the r.h.s. of ()3.5p 
and thanks to Assumption 12.31 we get for all stopping times 5^ < T„ < (5„ + 6) AT, that 



E(|^^f -^^f|)<5 ((p/(.,.)||oo + l)^ + ||/||oo(5 + d)J supE( sup 

^ ^ neN* te[o,T] 

+ ll/lloo^^ sup E( sup . (3.8) 

neN* te[0,T] 

The quadratic variation process (|3.6|) satisfies a similar inequality: 

E(|(M"-/)r„-(M"'/)5j) 



<II/IILe 



<ll/ll^5 



b{X-,l) , (i(X,",l) + t/(X,M) 



n n 



ds 



b+d\ ^, . ^, ,^2^ 



2f + — sup E( sup (Xf , 1)) + - sup E( sup {X"^ , ly) . (3.9) 
n ) nm* te[o,T] " neN* te[o,T] 



Then, for e > 0, > 0, a sufficiently large n and small (5, we have using ()2.1ip and Assumption 
that 

P(|A^f -A^f|>r?)<. and P(| (M"'/)t„ - (M"'/)5j > 7?) < e. (3.10) 

Prom (I22I), (I33I), (133]) and the Aldous-Rebolledo criterion (see e.g. [20] or [IH Th. 1.17]), we 
obtain the uniform tightness of the sequence ((X",/))ngN* in D([0, T],M). Thanks to Roelly's 
criterion [30j, we conclude that (X")„gN* is uniformly tight in B([0, T], (^M.p(X\v)). 

Let us denote by X a limiting process of [X'^^n&i* ■ It is almost surely (a.s.) continuous in 
{M.F{?i),v) since 

sup sup i(xr,/)-(xr,/)i<i. (3.11) 

/,ll/l|oo<l " 



In the case where is a compact subset of M'^, the vague and weak topologies coincide, which 
fails in the non-compact case. Nevertheless the tightness in B([0,T], {^M.p{X\w)) is needed to 
identify the limiting values of (X")„gN*. 

Step 2 Let us now concentrate on the case where X is unbounded and let us show the tightness 
of (X"')jigj^* in D([0,T], {^M.p{X),w)). The same computation as in Step 1 for /(x) = 1 implies 
that the sequence ((X", l))„g{!^* is uniformly tight in ID)([0, T], R+). As a consequence, it is 
possible to extract from {X'^^n&i* a subsequence (X"")„gN* such that: 

• (X"")„gN* converges in distribution to X in ]I])([0,r], {M.FiX),v)), 

• ((X"", l))ngN* converges in distribution in ©([O, T], M+). 

Let us now show that the limit of ((X"", l))„gN» is (X,l), empedding a loss of mass in the 
limit. Indeed, as a consequence, a criterion in Meleard and Roelly [25j will prove that (X"")„gN* 
converges in distribution to X in D([0,T], {^M.p(X\w)). 

By simplicity, we will again denote ti„ by n. 

As in Jourdain- Meleard ^21J, we introduce a sequence of smooth functions V'fc defined on and 
approximating l|„>;i.|. For /c S N, let tpk{u) = V'(0 V {u — {k — 1)) A 1) where tpi^y) = Qy^ — 15y^-|- 
lOy^ is a nondecreasing function such that -0(0) = ^"'(0) = '^"(O) = 1—^^(1) = ^'(1) = V'"(l) = 0- 
The function u i— )• ipk{u) is nondecreasing on M-|_, equals on [0, /c — 1] and 1 on the complement 



of [0,k). In particular ipQ = 1. Moreover the sequence {ipk)k£N* is nonincreasing, and satisfies 
for M > and p > 1 that 

^{u>k} < i'k{u) < l{„>fe-i}; (3.12) 
i;^^\u) < sup |V'i^)H|l|„>fc_i} < sup l^f (u)IV'fc-i(n). 

ue[k-l,k] ue[k~l,k] 

The proof of the following lemma is postponed at the end of Proposition I3.2[ s proof. We define 
fkix) = Tpk{\\x\\), for ah X £ X. 

Lemma 3.3. Under the assumptions of Proposition [XH 



lim limsupE sup{X^,fk) = 0. 

fc->+oo n-!>+oo \ t<T I 



From Lemma \'d.'6\ we can deduce that 

lim E(sup{Xt,fk)) =0. (3.13) 

Indeed, for A; G N, the continuous and compactly supported functions {fk/ *== /fe(l — fe))eeN 
increase to fk as i ^ +oo. By continuity of i— >■ supf^rpl^ut, fk/) on ]D){[0,T], {Mp,v)) and 
uniform integrability deduced from the uniform square moment estimates (|2.1ip . one has 

E(snp{Xt,fk,i)) = lim E(sup(Xr, A,,)) < Hminf E(sup(Xf , /fc)) . 

Taking the limit i — )• +oo in the left-hand-side by the monotone convergence theorem, one 
concludes that for A; = 

E(sup(Xt,l)) =E(snp{Xt,fo)) < +oo, (3.14) 
and from Lemma 13.31 that (|3.13|) holds for any k. 

As a consequence one may extract a subsequence of {sup^^j^{Xt, fk))k that converges to a.s., 
and since the process {Xt)t<T is continuous from [0, T] into (A1i;'(Af), t>), one deduces that it is 
also continuous from [0,T] into {AiF{'^)-,w). 

We can now prove the convergence of (X"",l) to (X,l). For F a Lipschitz continuous and 
bounded function on B([0, r],]R), we have 

lim sup |E(F((X", 1)) - F{{X, 1))) | < lim sup lim sup |E(F((X", 1)) - F((X", 1 - fk)))\ 

+ lim sup lim sup |E(F((X", 1 - fk)) - F{{X, 1 - fk)))\ 

fc— >+cxD n->+oo 

+ lim sup |E(F((X, 1 - fk)) - F{{X, 1))) |. 

k—^+oo 

Since |-F((j^, 1 — fk)) — F{u, 1))| < C sup^<'p(ft, /fc) by Lipschitz property, the first and the third 
terms in the r.h.s. are equal to respectively according to Lemma [3.3l and to ()3.13p . The second 
term is by continuity oi v ^ {u,l - fk) in ]I])([0,r], {Mf{X),v)). 

This ends the proof of Proposition 13.21 ■ 
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Proof of Lemma \3.3l Firstly, let us show that for each t G [0, T], 

lim limsupE = 0. (3.15) 

The boundedness of r and Assumption 12.31 -2 ensure the existence of a sequence (en)raeN* con- 
verging to such that 



Jo Jo 
+ e( / / r{x,a)p{x,a) Afk{x)X'^{dx,da)ds) , (3.16) 

and we have by and (f3J2]l 

/ r{x,a)p{x,a)Afk{x)Xl\dx,da) < f \\D'fk\\^ {X^Jk-i). 

Since moreover, the sequence {fk)k<^N* is non-increasing, < {X^,fk-i) and there is a 

constant C > independent of A; > 2 such that 

E((X^/fc)) <E((Xo",/fc)) + C fE{{X^Jk-i))ds + en f E{{X^,l))ds. (3.17) 

JO JO 

Let /i"-'' = E ((X,", /fc)) < /u^ = E ((X,", 1)) which is bounded uniformly in n e N* and s E [0, T] 
according to (|2.1ip . There exist two positive constants Ci and C2 such that 



rt 

n,k ^ n,k , ^ / n.fc— 1 j , /-i 

lit < ^0 +'^1 / ds + C2en. 

Jo 



Iteration of this inequality yields 



£=0 ■ ■ £=0 

£=Lfc/2j+i ■ ^ '' 

where we used the monotonicity of /Iq''^ w.r.t. k for the second inequality. Given the moment 
condition (|2.1U|) . the assumption of tightness in (7Wir(A'),w) of the initial conditions (Xq)^^^* 
is equivalent to 



n,k 

illli lllli SUJJ /-' 

fc->+oo n_>+oo 

Hence 



lim limsup/UQ* =0. (3.18) 

fc->+oo n_>+oo 



iim lim sup //j < sup //q hm > — h iim —TT-rr- ■ 

k^+oo n^+oo nSN* A;-5>+oo f-^ V. k-^+oo [k)\ 

l=\k/2\+l 



We deduce immediately that 



lim limsupE ((Xf = lim limsupAi^ =0. (3.19) 

fc-s>+oo „_!.+oo A:-5-+oo n-5-+oo 

Let us now consider the martingale M^'^ defined by ()3.5p with fk instead of /, and with 
quadratic variation given in ()3.6p . Similar arguments as above allow us to prove that 

E((M"''=)t) < Ci [\{{X^Jk^i))ds + enC2 fE{{X'^,l))ds. 
Jo Jo 
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Thus, using that < 1, Doob's inequahty, (|3.15p . (j2.1ip and the dommated convergence 
theorem, we get 

hm hm sup E ( sup | M"''^ | ) = 0. 

k^+oo n-s>+oo ^ t<T ' 

Let us now come back to the process (^", /fc)- As before, we can get 

(Xr, h) < {X^, fk) + Ml"'" + b [\x2, fk)ds + En ^)ds 

Jo Jo 

+ / r{x,a)p{x,a) Afk{x)Xg{dx,da)ds 

Jo JR'ixR+ 

< {XS, fk) + M^''' + Ci [\x^, h-i)ds + En C2 Ax,", l)ds, (3.20) 

Jo Jo 

for constants Ci and C2. Let a"'^ = E (sup,<((X", /fc)) and a" = E (sup,<((X", 1)) which is 
bounded uniformly in n S N* and t G [0, T] according to ()2.1ip . One deduces that 

^n,k ^ ^n,k ^ /■* ^n,fe-l^^ ^ ^^^^ + SUp IM,"''^!) . 

An iteration as before ahows us to prove that 

hm hmsupE( sup(X", /fc) ) = hm hmsup a"''^ = 0, 



fc-s>+oo n-s>+oo ^t<T ' fc->-+oo n-!>+oo 

which concludes the proof of Lemma 13.31 and thus the one of Proposition [37 



3.2 Identification of the limiting values 

To obtain the convergence stated in Theorem 13. H we show that the limiting value X of the 
uniformly tight sequence (X")„gN* is unique. We establish a martingale problem satisfied by 
X in which there are integration terms with respect to the equilibrium ()3.2p involved in the 
averaging phenomenon for the ages. The uniqueness of the solution to the martingale problem 
is then proved. 



3.2.1 Averaging phenomenon 

We begin with establishing the form of the limiting values of the time-marginal distributions 
(X"((ix, da))„gN* for t G [0,T]. Since the sequence (X"')„gN. is uniformly tight, there exists a 
subsequence of (X"((ix, (ia))„gN* , with trait-marginals converging in law to a limiting value X, 
that by simplicity, we denote again by {X'^{dx^da))n&\* ■ 

We have already explained why the uniform tightness of the sequence {t ^ X^idx^ da))n£N* in 
O{[0,T],AiF{X X M_|_)) cannot hold. However, following Kurtz [23], we will prove the uniform 
tightness of the sequence of random measures 

r"((it, dx, da) = X^{dx, da)dt (3.21) 

on A^_f([0, T] X A" X M+). Proceeding in this way allows us to escape the difficulties created by 
the degeneracies due to the rapid time scale for age, when one tries to follow individual paths. 

Proposition 3.4. Under the assumptions of Theorem \3.1\ the sequence (r")„gN* converges in 
law to Xt{dx)m{x,da)dt in Mf{[0,T] x ^ x R+). 

As a consequence, for dt- almost every (a.e.) t £ [0,T], the sequence {X^{dx, da))n£n* converges 
weakly to rh{x,da)Xt{dx), with fh defined in i3.2\) . 
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The proof of Proposition 13.41 is inspired by Kurtz [23]. To establish the the result, we need to 
consider the random measure defined in (|3.2ip . 

Firstly, we prove the uniform tightness of the sequence {Xl^{dx,da))nef'i* , ^or fixed t £ [0,T] 
(Lemma 13. 5p . as well as the one of the sequence of measures (r")„gp^* fLemma l3.6p . where the 
pathwise and individual points of view have been forgotten. The techniques to disentangle 
the traits and individuals' time scales appear strikingly in the proof of Lemma \3.5\ where 
different treatments are used for the trait marginal and for the ages, with the introduction 
of the individuals' lifelengths. Then, in the proof of Proposition 13. 4t a factor n appears in 
()3.39p . when changing from the macroscopic scale to the microscopic scale. The next part of the 
proof consists in identifying the limiting martingale problem. 

Lemma 3.5. For dt-a.e. t G [0,T], the sequence {X^)nen* is uniformly tight on Aip{X x M-|_). 

Proof. Let e > 0. Since the family (X")„gp|* is tight, there exists a compact set K C M*^ such 
that 

sup ¥{XJ'{K^) >e) <e. (3.22) 

neN* 

Moreover, thanks to Point 2 of Assumption 12.11 and using a coupling argument, the life-lengths 
of the individuals in the population X" born after time are dominated, uniformly in x £ X, 
by independent random variables Df (t) with survival function 5" defined in (|2.4p . Because the 
aging velocity is n, the ages of these individuals satisfy Ai{t) < nD^{t). For an individual i 
alive at time 0, conditionally on the state at time 0, the remaining time it has to live can be 
dominated by an independent random variable A" such that 

P(A^ > e) =exp(^- nr{Ai{0) + nu)du^ . 
Thus, for ^ > and no, iV G N*, 

1 ^" 

sup F{X^i{K X [0,^])'=) > 3e) < sup F{X^{K'^) > e) + sup p(- J] l{A,(t)>A} > 2e) 



n>no n>no n>no i—1 

nN 

<e + sup ^{«^r(t)>^} > + ^^^P ^(^" > 

n>no ^ ^ n>no 

+ sup P( ^ '^{Aao)+nA^>A} > nej . (3.23) 

n>no ^ -^^ ' 

By (j2.10p . it is possible to find such that: 

sup P(iVr > nX) = sup P((Xr, 1) > iV) < ^^Pn>noE(^UP..[O.T](^r,l)) ^ ^ ^3 ^4) 

n>no n>no 

Let us consider the second term of (|3.23p . Notice that 

IE(l{„D^-(t)>A}) = S'^^Ajn) = exp - / r{a)daj , 



which converges to when A tends to infinity by ()2.3p . For as in (|3.24p . there hence exists 



13 



Ai sufficiently large so that for A> Ai, S'^{A/n) < e/2N. Then 

nN nN 

^(E ^{^rW>^M > = IP(E(%rW>^M - ^''{A/n)) > n{E - iV5"(^/n)) 

i=l i=l 

< P( J^(]l|B^n(,)>^/„|-S'^(A/n)) > ne/2) < exp 



i=l 



(3.25) 



by Bernstein's inequality {e.g. [32] p. 855). For a sufficiently large n\ and for n > rii, the r.h.s. 
of ()3.25p is smaller than e. 

Let us now upper bound the last term of (|3.23p . Let A 1— )• (piA) be such that (p{A) and 
Ii^{A) L{u)du tend to infinity when A tends to infinity. We firstly notice that: 

<Xo"(Af X [</.(A), +00)) + (Xo", l)e- (3.26) 

Recall that Xq is the limit, in probability and for the weak convergence, of (Xq)^^^*- There 
hence exists n2{A) sufficiently large such that for n > n2{A), the right hand side of (|3.26|) is 
smaller than 

X^{X X [<A(A),+oo)) + (Xo,l)e~^^1^)^(")''" + I 

with probability 1 — e. By choice of (t){A), there exists A2 sufficiently large such that such that 
for A > A2, and n > 722(^2) this upper bound is smaller than e/2. For such A and n, we have 

1=1 i=l 

2 

<E f exp f - -) )+e. (3.27) 

by applying Bernstein's inequality again. There exists 723 sufficiently large so that for every 
n > n^, the right hand side of ()3.27p is smaller than 2e. 

The tightness of (Xj")„eN* is thus a consequence of IK2^ . IKMi . IK2^ and ([3:271) . with the 
choices of j4 > max{Ai, A2) and no > max(ni, 712, ^3). ■ 

Lemma 3.6. The family (r")„eN is tight inMpiiO^T] x X x M+). 

Proof. Following Kurtz |23[ Lemma 1.3], a sufficient condition for the tightness of the family 
(r")„gN is that for all e > 0, there exists a compact set H of x IR+ such that 

sup E(r'^([0,r] X E^)) < C{T)e. (3.28) 

Let us establish ([3.28p . From the proof of Lemma (331 it appears that the upperbounds ([3.22p . 
([3.23P and (|3.25p are uniform in t € [0, T] so that: 

sup sup f(x^{{K X [0,A]y) > 3e) < 5e. (3.29) 
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We are now ready to upperbound 



E(r"([0,r] x{Kx %A]r)) =e(^ (Xr,%x[o,A])=>di) = ^ IE(^r((i^ X %A]r))dt. 
Indeed: 

E(xf ((i^ X [0,^])'=)) < 3e Pfe((i^ X [0,^]r) < 3e) + E((Xr, I)! 



•'[(Kx[Q,A\Y)>^e 

<3e + yE((Xr,l)2)yp(x,"((i^ X [0,^])^=) >3e) <C7(r)(e + V^), (3.30) 
by Cauchy-Schwarz inequality and (|2.1ip . This proves p.28p and finishes the proof. ■ 

Before proving Proposition 13.41 we provide a lemma characterizing m(x,a). 

Lemma 3.7. Let x £ X be fixed. There exists a unique probability measure fh{x,da) on 
solution of the following equation: For ip € C^(M+,M) with compact support in [0,+oo), 



daip{a)fh{x,da) = / 'ip{a)r{x,a)fh{x,da) — ip{0) / r{x,a)fh{x,da). (3.31) 

The probability measure fh{x, da) is absolutely continuous with respect to the Lebesgue measure 
and its density is given in i3.^} . 

Proof. Let us consider the test function ip{a) = J^°° f{a)da, where / G Cc(lR+,I^+) is positive. 
Then daip{a) = —f{a) and V'(O) = /o^°° f{a)da. Equation ()3.3ip gives by Fubini's theorem: 



f {a)rh{x , da) = j J f{a)dar{x,a)fh{x,da) 

r+oo 

f{a) / r{x,a)m{x,da) da. (3.32) 



This entails that m(x, da) is absolutely continuous with respect to the Lebesgue measure with 
density m{x,a) = J^°° r{x,a)fh{x,a)da. The latter implies that a i— )• m(x,a) is a function of 
class C^. Using further an integration by part in ()3.3ip . and since m(x,a) tends to when a 
grows to infinity, we get for all ip £ Cl (1R+ , M) 

-mrn(.. 0) - f = [ (Ma) - m) r(^, a)B(.. a)d„. (3.33) 

By identification, we obtain that m(x, a) is a solution of 

dafh{x, a) = —r{x, a)fh{x, a) 

fh{x,0) = / r{x,a)fh{x,a)da, (3.34) 

which is solved by 

m(x, a) = m(x, 0) exp ( — ^ r{x,a)da). (3.35) 
Since fh(x, a)da is a probability measure, necessarily 







m{x,0) = r{x,a)fh{x,a)da = r^^— . , ^ , . (3.36) 

Jr+ Ju+ ( - Jo r{x,a)da)da 

This provides existence and uniqueness of the solution of (j3.34p and hence of ()3.3ip . ■ 
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Remark 3.8. Notice that the system (j3.34p defines the stable age equihbrium of the McKendrick- 
Von Foerster equation [24j il7j (see also [36j) when the birth and death rates equal to r{x,a) 
and the trait x is fixed. □ 



Remark 3.9. The space C^(M+,M) is separable and there exists a denumberable dense family 
{ipk)k<^N in this set. To obtain the result of Lemma |3.7| it is sufficient to have (|3.31|) for these 
functions {'ipk)keN- 



We are now able to prove Proposition 13? 



Proof of Proposition [X^ From (|3.2ip . we can see that the marginal measure of r'^{ds,dx,da) 
on [0,T] X is X'^{dx)ds. For any real bounded test function ip : (s,x) i— >• ips{x) on [0,T] x Af, 



/7 

Jo J A 



ipsix)T"-{ds,dx,da) = / {X^,Lps)ds. 



(3.37) 

Ja:xr+ Jo 

The sequence {X"')nm* is uniformly tight by Proposition [321 as well as (r"')„gN*, by Lemma [3^ 
(ii). Using Prohorov's theorem, we thus deduce that {T^{ds,dx,da),X'^{dx)ds)n is relatively 
compact and there exists a subsequence that converges in distribution to a limiting value, say 
{T{ds,dx,da),Xsidx)ds). Taking ()3.37p to the limit, we obtain that Xs{dx)ds is necessarily 
the marginal measure of T{ds, dx, da) on [0, T] x up to a null-measure set. We deduce from 
this {e.g. Lemma 1.4 of Kurtz [23]) that there exists a (random) probability-valued process 
{7s,x{da), s £ [0,T],x G X) that is predictable in (a;, s) and such that for all bounded measurable 
function if{s,x,a) on [0,T] x X x M_|_, 



if{s, x, a)T{d.s, dx, da) 



JXx 



JX 



ip{.s, X, a)^s,x[da)Xs[dx)ds. 



(3.38) 



We now want to characterize the limiting value T(ds,dx,da) = ^s,x{da)Xs{dx)ds. Applying 
(|2.12p for a test function <p{x,a) G C^'^{X x 



n 



{x^,^)-{xs,v) 



n 



+ pix, a) (^r{x, a) + 



JXxR+ 

b{x, a 



n 



b(x,a) , , ( d(x , a) + X'}\J (x , a 



n 



^) and dividing by n gives that: 

daf{x, a) + r(x, a) ^v^(x, 0) — '^{x, a) 
(yip{x + /i, 0) - (^(x,0))7r"(x,a,d/i) 

(/?(x,a)lr"(ds,dx,da) (3.39) 



n 



is a martingale. Using ()2.1ip . we can easily prove that (M"''^/n)^gpj, te[or] i^ uniformly inte- 
grable and that 



lim E 

n— >+oo 



n 



where 



Mf 



JXx 



da^{x, a) + r{x, a) [ip{x, 0) — ip{x, a)) js,x{da)Xs{dx)ds 



(3.40) 



Moreover, the uniform integrability of (M'^'^/n) also provides that the process {Mf)t is a 
martingale. As it is also a continuous and finite variation process, it must hence be almost 
surely zero. Since this holds for every t E M+, we have proved that a.s., dt-a.e. 



/ da(p{x,a) +r{x,a){ip{x,0) - (p{x,a)) jt,x{da)Xt{dx) = 0. 

JxxR4- '- 



(3.41) 
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Choosing ip{x, a) = (l)i{x)il^k[a) with {4>i)(,^n and {il^k)k&i dense famihes in Cc{X , M) and C].{ 
respectively, we obtain from (I3.4ip that, for all G N, 



c^i{x)Hk{t,x)Xt{dx) = {) 



X 



where Hk{t,x) = J^^ daipkio-) + i"{x,a)[ipk{0) — V'fc(fl)) lt,x{da). Almost surely, the function 

Hk{t,x) is bounded and is thus dt-a.e. Xt((ia;)-integrable. We obtain that for all A; G N, a.s., 
dt-a.e., Xt{dx)-a.e., 



datpk{a) + r{x,a){^pk{0) - ipkia)) 7t,xida) = 0. 



(3.42) 



By Lemma 13.71 and Remark 13.91 we deduce that a.s., dt-a.e., and Xtidx)-a.e., ^t,x{da) = 
fh{x,a)da and as a consequence, any limiting value of (X")„gN. is of the form Xt{dx) (gi 
m{x,a)da. ■ 

3.2.2 Characterization of the limiting values 

In the previous sections, we have proved that the sequence (X")„gN* is tight and that for a given 
limiting value X, the associated subsequence {T^{dt,dx,da) = X^{dx,da) dt)nm* converges in 
(A^i?([0,T] X A' X R4.),?i)) to Xt{dx)m{x,da) dt. Now, we are ready to prove that: 

Lemma 3.10. The limiting values X of the sequence (X")„gp^* are solution of the martingale 
problem (Egj-fOj. 



Proof Let < si < . . . Sk < s < t, and let us introduce for Y G D(M+, A^f('^)): 



,,t{Y) = MYs,)...MYs,){{Yt,f) - {Ys, f) - I' du I^Yuidx) 



(^){x)Af{x) + (b{x) - d{x) - YsU{x))f{x)\ }, (3.43) 



where (/>!, . . . ,(f)k are bounded continuous functions on Aip(X) and / € T^{A). Our purpose is 



to prove that E(^s^f(X)) = for any limiting value X of {X^] 



neN* 



Let X be a limiting value of {X^)ni=N* and let (X"")„gp}* be a subsequence converging to X. 
On the one hand, thanks to Proposition 13.41 ()2.1ip and (|3.1ip : 



E{^sAX))= lim E(MK")---MK"){(^t"J)-{K"J)- fdu [ 



X^"{dx,da) 



p{x,a)r{x,a)Af{x) + (6(x, a) - d(x, a) - / Ui{x, a), {y,a))X:^" {dy,da)) fix)]}) . 

(3.44) 



On the other hand, the term under the expectation in the r.h.s. of (|3.44p equals: 



(3.45) 
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where M""'-^ has been defined in (13.51) and where: 



du 



t 

du 



X^"{dx,da) r{x,a) 



X 



fix + h)- f{x)j iv:"" (x, a, dh) - p{x, a)Af{x) 
f{x + h)-f{x)y^"{x,dh)-Afix) 



du 



X^"{dx,da) r{x,a)p{x,a) 

Xl^ (dx, da) b{x, a) \ [ f{x + h)K''" (x, a, dh) - f{x) 



Firstly, using ()2.1ip and the fact that the process M"''^ is a martingale we obtain that: 

E(0i(X,7)...,^,(X«;)[M,"-^-M^'^]) =0. (3.46) 

Secondly, from Assumption 12.31 \ f{x + h)K'"-{x,a,dh) — f{x)\ = o{\/n) and using (|2.11|) 
again provides: 

liin E(</.i(X,7)...0fc(X,";)[A,„+i?J) =0. (3.47) 

From ([333]), dOSjl . ^Mh and (f3:i7D . we deduce that E{^s,t{^)) = and hence {Mf)tm+ 
defined in ()3.3p is a martingale. From this, using Ito's formula with localization arguments and 
Proposition 12.61 (i). we obtain that the following process is a martingale: 



{Xt,f)''-{Xo,f)' 



2{Xs,f) I {pr){x)Af{x) + (bix) 

- d{x) - XsU{x))f{x)]Xs{dx) ds - {M^)t (3.48) 



X 



is a martingale. Moreover, using the results of Proposition \2.Q\ we obtain that the following 
process is a martingale: 



(xr,/)2-(Xo",/)2 







Xxl 



nr{x,a) {f{x + h)-f{x))K''(x,a,dh) 



JXx 



+ b{x,a) / f{x + h)K''{x,a,dh)-{d{x,a)+X'^U{x,a))f{x))Xl'{dx,da) 
r{x,a)( [ f{x + h)K''{x,a,dh) + f{x)) 
f{x + h)K''{x,a,dh) + 



ds 



+ 



b{x, a 



\x + h)K-{x^ a, dh) + di^^^)+XsUix , a) 
n J^d" n 



f\x) X^{dx,da) ds. 



By using arguments similar as those in the beginning of the proof, we deduce that 

rt 



2{Xs,f) / {pr)Af{x) 

'X L 



{Xt,f)^-{Xo,ff- 

+ (b{x) - d{x) - XsU{x))f{x) Xs{dx)ds 



2?{x)f{x)Xs{dx) ds (3.49) 



JX 



is a martingale. Comparing (|3.48p and (|3.49p yields the bracket of the martingale . This 
ends the proof. ■ 
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3.2.3 Uniqueness of the martingale problem 

We have shown that the hmiting values of the uniformly tight sequence (X'^)n^^* satisfy the 
martingale problem (|3.3p - (|3.4p . To conclude the proof of Theorem 13.11 it remains to prove the 
uniqueness of the solution of this martingale problem. 

Proposition 3.11. There is a unique solution to the martingale problem of Theorem \3.1\ 

Proof. We start with getting rid of the non-linearity by using Girsanov's formula (see Dawson 
[8j Theorem 7.2.2). There exists a probability measure Q on the path space such that for all 
feV{A): 

Ml =mI + [ [ (Kx) - d{x) - XsU{x)) f{x)Xs{dx) ds 

Jo jR'i ^ ' 

= {XtJ)-{XoJ)- f [ (P){x)Af{x)Xs{dx)ds (3.50) 
Jo JR'^ 

is a square integrable martingale with bracket ()3.4p . 

The uniqueness of the solution the martingale problem (|3.50p - (|3.4p is proved by Roelly and 
Rouault [3lj . It is based on the branching property of X under Q which allows us to characterize 
the Laplace functional Ltf of Xt by its cumulant Utf: 

Li(/) = E(e<^'/)) =E(e<^0'^*^)). (3.51) 

The latter is the unique positive solution of the following PDE: 

— {t,x) = Au{t,x)-?{xy{t,x), u{0,x) = f{x), (3.52) 

(see e.g. Pazy [Ml Th. 1.4 and 1.5 p.l85 and 187]). 

From the uniqueness of the solution of ()3.50p - ()3.4p . we deduce classically the uniqueness of the 
solution of (|3.3p - (j3.4p . (See for example Evans-Perkins |12j, Fitzsimmons p[B])- B 

The proof of Theorem 13.11 is now complete. 

4 Examples 

Let us develop and compare two examples, which only differ by the function r{x, a). 
4.1 Logistic physical-age and size-structured population 

In Meleard and Tran [26| , the following example for a population structured by age a G M-(_ and 
size X G X = [0, xq] is considered: 

6(x,a) = x{xo - x)e""l[o,2^(,](x) for xq > 0, 

d{x,a)=do, U{{x,a),{y,a)) = r]{xo - x), (4.1) 

with Xq = 4, do = 1/4 and r] = 1.7. Because reproduction needs energy, and since this energy 
depends on the size of the created offspring, very small or big individuals are disadvantaged. 
Individuals of intermediate size x = 2 have the highest birth rate. The competition term 
in contrast favors bigger individuals. Hence there is a trade-off between competitiveness and 
reproduction. The decreasing exponential in age describes a senescence phenomenon: older 
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individuals reproduce less than their young competitors. In [26], partial differential equation 
limits, Trait substitution sequence and Canonical equations are considered. Here we consider 
the superprocess approximation described in the above sections, with r{x,a) = 1 and ■7T"'{x,dh) 

2 

a centered Gaussian kernel with variance ^ conditioned on [0, xq], as in Example 12.51 

Computation gives fh{x,a) = e~°' so that Xt{dx,da) = Xt{dx) ® e~"' da becomes in this 
particular case a product measure. As soon as the population survives, the age distribution 
"stabilizes" around an exponential distribution with parameter 1, as seen on the simulations of 
Figure [TJ With the age distribution m{x,a) = e"'*, we get 

b{x) = x{xo -x) [ e-^^da = "^^^^ ~ ; d{x) = do U{x, y) = r/(xo - x). (4.2) 



The martingale problem f|3.3p becomes here: 

m/ = {Xt, f) - (Xo, f)- (pyA/(x) 

+ [^^^ - (do + vixo - x){Xs, l))]f{x))Xs{dx) ds, (4.3) 

{Mf)t= f I 2f\x)Xs{dx)ds. 
Jo Jx 

In Figure [U two sets of simulations are presented, depending on two different mutation variances 
cr^. As expected, when a increases, the traits vary more rapidly, and the irregularity of the 
trait support appears more strikingly. On both simulations of Fig. [H extinction happens in a 
fast time. Almost-sure extinction is due to the logistic interaction, as proved in the following 
proposition. 

Proposition 4.1. There is almost-sure extinction of the superprocess ( [^.3[ ). 
Proof. The mass of the super-process satisfies the following equation: 

(Xt,l) = (Xo,l)+ f [ A{x,{Xs,l))Xs{dx)ds + M} 

Jo J[0,xo] 



{M%= [\{XsA)ds, 
Jo 



where A{x, Z) = "^^^"^ - {do + r]{xo -x)Z). (4.4) 

This equation is not closed for the mass process, since the drift depends on the trait distribution. 
Our purpose is to upper-bound A(x, Z) so that (X, 1) can be stochastically dominated by a Feller 
diffusion with negative drift, that goes extinct almost surely. 

In the case where x > xq — (^ — C) with C, £ (0, ^ A 1) and since x £ (0, xo), one gets 

A(x, Z) = - do + {xo ~ ^) (I ~ 

<-do + (— -C)x- = -— . (4.5) 



In the case where x < xq — (^ — C), then < ^ — C ^ xo — x < xq and depending on the sign 
of x/2 -r]Z: 



A{x,Z) < -(io + max(^xo(y - vZ) ; ~ ( y ~ 

<f-do-v{—-C)Z. (4.6) 

/ Xo 
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Figure 1: Simulation of the individual-based process X", with n = 1000 and discretization step 
At = 0.005. The system is started with 1000 particles of trait x = 1.5. First line: a = 1. Second 
line: a = 0.8. (a): Support of the process X" (with time in abscissa and trait in ordinate), (b): 
Evolution of the population size, (c): Age distribution for t = 0.5. It can be checked that the 
age distribution converges to an exponential of parameter 1 (plain line). 



"^0 = ' (4.7) 



Since the upper bounds in ()4.5p and (j4.6p are equal when the mass Z equals rriQ defined by 

we thus get in any case that 

A(x, Z) < -^lz>mo + (4 -do) lz<mo- (4.8) 



2 - V 2 
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Hence, the process (X , 1) can be stochastically dominated by the following positive process: 



Zt = {Xo, 1) + ^ ( - + ^p( ^°(^o + _ do)lz.<m,)ds + ^ V^sdB, (4.9) 

where i? is a standard Brownian motion. 

We can adapt the results of Meyn and Tweedie j27j to prove almost sure extinction. For 
u £ M+, let us denote by = inf{t > 0, Zt < u} and let z = {Xq,1) > 0. Either z < mo 
and then Tmo = 0, or z > uiq. In the latter case, let us consider M > z > mo and let 
= mi{t > 0, Zt> M}. We have 

Y 



( ^r„, APM - ^ + V / Zsds- I yYz.dBs ) = 0. (4.10) 



Jo 



By the uniform integrability of the fourth term and the optional stopping theorem, we have 

Since moreover Zg > mo for s G [0, A pu], we deduce that 



aTruQ^pM V 2Z 

Zs ds) < —. (4.11) 

It can easily be proved that for all T > 0, E(supj<y(Zj)^) < oo, implying that pM tends to 
infinity with M. Thus, ()4.1ip provides that for all z > 0, Pz(rmo < +oo) = 1. By Girsanov's 
theorem, there exists a probability measure under which the process Z is a sub-critical Feller 
diffusion. It turns out that Pmo(''"o A pM < +oo) = 1. Standard computation using the strong 
Markov property yields P^(to < +oo) = 1. ■ 

4.2 Logistic biological-age and size-structured population 

In this section, the trait x £ [xi, X2] C (0, xq) (with xi,X2 > 0) is linked to the rate of metabolism, 
which measures the energy expended by individuals, and is often an increasing function of the 
body size. Ageing may result from toxic by-products of the metabolism. This leads us to 
introduce a biological age xa, where a is the physical age and where x can be interpreted 
as the ageing velocity equals. In this example, we consider r(x,a) = xa so that biologically 
older individuals give birth and die with higher rate, the other parameters being chosen as in 
Subsection 4.1. For a review on body size, energy metabolism and ageing, we refer the reader 
to [33]. In our example: 

xa'^ 

m[^x,u) = y= — ]l[o,+oo)(o)- (4.12) 

V 27r 



We recognize the Gaussian distribution with variance 1/x conditioned on being positive. Then: 

, 2x3/2(xo-x) /■+°° .(a2+2a/.) 2x(xo-x)e'^ /■+°° , 

b[x) = -= / e 2 da = — ^ — / e 2 da 

V27r Jo v27r Ji/y^ 

=2x{xo-x)e^^( --^), (4.13) 



^ , . 2:e3/2 /■+- _^ /2^ /■+°° 

r[x) = — / ae ^ da = \ — / ae ^ da = \ — , (4-14) 




2tt Jo V vr JO 

where $ is the cumulative distribution function of the standard Gaussian distribution. The 
functions d{x) and U{x,y) are unchanged and given by ()4.2p . The martingale problem ()3.3p 
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becomes here: 




(4.15) 



(b) 



Time 1.500001 Time 0.5000002 




Figure 2: Simulation of the individual-based process X"^ , with n = 1000 and discretization step 
At = 0.005. The system is started with 1000 particles of trait x = 1.5. First line: a = 1, Second 
line: a = 0.2. (a): Support of the process X" (with time in abscissa and trait in ordinate), (b): 
Evolution of the population size, (c): Marginal age distribution for t = 0.5. For comparison, we 
draw the density iri{l, a) (plain line), m(0.5, a) (dotted line) and iri{3,a) (dashed line). 



In this example, there is a higher senescence for individuals with trait x > 1, compared with the 
example of Section [4.11 The new choice of r(x, a) influences the age distribution: lifelengths are 
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Figure 3: Joint distributions in trait and age for the simulation of Fig. [21 We see that contrarily 
to the example of Section \4.1[ we do not have here independence of traits and age. (a): a = 1. 
(h): a = 0.2 



shortened. This can be seen on the smaller support of the age distribution (compare Fig. [2]-(c) 
with Fig. [I]-(c)). 

However, the populations are more persistent in the example of this section, although it can 
be proved similarly to Prop. 14.11 that there is almost sure extinction. Indeed, contrary to the 
populations in Example 1 which are extinct at t = 2, the population of Example 2 still survives 
at t = 20. One reason is that the growth rate in the finite variation term of (|4.15|) is bigger than 
the one in Example 1. Indeed, for many values of the factor 2exp(l/2x)<I>(— l/y^) in the 
birth rate h{x) is bigger than the factor 1/2. For x = 1.5, 2 exp(l/2x)<I>(— 1/^/x) = 0.58 > 0.5 
and for x = 3, 2 exp(l/2x)^>(-l/i/x) = 0.67 > 0.5. 

When comparing Fig. [I]-(b) and Fig. [2]-(b), we observe more fluctuations of the population 
size in Example 2. The bracket of the martingale in ()4.15p presents a multiplicative x term, 
compared to (|4.3p . As soon as x > |, this explains the increased variance. Notice however that 
this variance tends to zero when the population size tends to zero, which also explains why there 
is no decrease in the population persistence. 

Finally, the multiplicative term p^x/2-Ka'^ in front of the diffusion term A(/p(x) explains the 
large variability of the trait support, which is observed in Fig. [2]-(a). When the diffusion coef- 
ficient a is small (second line of Fig. [2D, the traits evolve towards a value between x = 2 and 
X = 4 where the trade-off between reproduction and competition is optimized. 

Acknowledgements: This work benefited from the support of the ANR MANEGE (ANR- 
09-BLAN-0215) and from the "Chaire Modelisation Mathematique et Biodiversite of Veolia 
Environnement-Ecole Polytechnique- Museum National d'Histoire Naturelle-Fondation X". The 
authors thank S. Billiard, T. Kurtz and L. Popovic for enriching discussions. 

References 

[1] K. Athreya, S. Athreya, and S. Iyer. Super-critical age dependent branching Markov processes and their 
scahng hmits. Bernoulli , 17(1):138-154, 2011. 

[2] K. Ball, T.G. Kurtz, L. Popovic, and G. Rempala. Asymptotic analysis of multiscale approximations to 
reaction networks. Annals of Applied Probability , 16(4):1925-1961, 2006. 

[3] A. Bose and I. Kaj. Diffusion approximation for an age-structured population. Annals of Applied 
Probabilities , 5(1):140-157, 1995. 



24 



[4] A. Bose and I. Kaj. A scaling limit process for the age-reproduction structure in a Markov population. 
Markov Processes and Related Fields , 6(3):397-428, 2000. 

[5] N. Champagnat, R. Ferriere, and S. Meleard. Unifying evolutionary dynamics: from individual stochastic 
processes to macroscopic models via timcscalc separation. Theoretical Population Biology , 69:297-321, 2006. 

[6] N. Champagnat, R. Ferriere, and S. Meleard. From individual stochastic processes to macroscopic models 
in adaptive evolution. Stochastic Models , 24:2-44, 2008. 

[7] B. Charlesworth. Evolution in Age structured Population . Cambridge University Press, 2 edition, 1994. 

[8] D. A. Dawson. Mcsurc-valucd markov processes. In Springer, editor, Ecole d'Etc de probabilites de 
Saint-Flour XXI , volume 1541 of Lectures Notes in Math. , pages 1-260, New York, 1993. 

[9] D.A. Dawson, L.G. Gorostiza, and Z. Li. Nonlocal branching superprocesses and some related models. Acta 
Applicandac Mathematicae , 74:93-112, 2002. 

[10] E.B. Dynkin. Branching particle systems and superprocesses. Annals of Probability , 19:1157 1194, 1991. 

[11] A. Etheridge. An introduction to superprocesses , volume 20 of University Lecture Series . Providence, 
American Mathematical Society edition, 2000. 

[12] S.N. Evans and E.A. Perkins. Measure-valued branching diffusions with singular interactions. Canadian 

Journal of Mathematics , 46:120-168, 1994. 

[13] S.N. Evans and D. Steinsaltz. Damage segregation at fissioning may increase growth rates: A superprocess 
model. Theoretical Population Biology , 71:473-490, 2007. 

[14] R. Ferriere and V.C. Tran. Stochastic and deterministic models for age-structured populations with geneti- 
cally variable traits. ESAIM: Proceedings , 27:289-310, 2009. 

[15] P.J. Fitzsimmons. On the martingale problem for measure-valued markov branching processes. In Seminar 
on Stochastic Processes (Los Angeles, CA, 1991) , volume 29 of Progr. Probab. , pages 39-51, Boston, 1992. 
Birkhaiiser. 

[16] K. Fleischmann, V.A. Vatutin, and A. Wakolbinger. Branching systems with long-living particles at the 
critical dimension. Theoretical Probability and its Applications , 47(3):429-454, 2002. 

[17] H. Von Foerster. Some remarks on changing populations. In Grune & Stratton, editor. The Kinetics of 
Cellular Proliferation , pages 382-407, New York 1959. 

[18] P. Jagers and F. Klebaner. Population-size-dependent and age-dependent branching processes. Stochastic 

Processes and their Applications , 87:235-254, 2000. 

[19] P. Jagers and F. Klebancn. Pojjiilatioii size (iepciidcnt, age structured branching processes linger around 
their carrying capacity. Joiuiial of Applit'd Probability , page to appear, 2011. 

[20] A. Joffe and M. Metivier. Weak convergence of sequences of semimartingales with applications to multitype 
branching processes. Advances in Applied Probability , 18:20-65, 1986. 

[21] B. Jourdain, S. Meleard, and W. Woyczynski. Levy flights in ecology, preprint HAL-00560633. 

[22] I. Kaj and S. Sagitov. Limit processes for age-dependent branching particle systems. Journal of Theoretical 
Probability , ll(l):225-257, 1998. 

[23] T.G. Kurtz. Averaging for martingale problems and stochastic approximation. In Springer, editor. Applied 
stochastic analysis (New Brunswick, NJ, 1991) , volume 177 of Lectures Notes in Control and Inform. Sci. , 
pages 186-209, BerUn, 1992. 

[24] A.G. McKendrick. Applications of mathematics to medical problems. Proc. Edin. Math.Soc , 54:98-130, 
1926. 

[25] S. Meleard and S. Roelly. Sur les convergences etroite ou vague de processus a valeurs mesures. 
C.R.Acad. Sci. Paris, Seric I , 317:785 788, 1993. 

[26] S. Meleard and V.C. Tran. Trait substitution sequence process and canonical equation for age-structured 
populations. Journal of Mathematical Biology , 58(6):881-921, 2009. 

[27] S.P. Meyn and R.L. Tweedie. Stability of Markovian processes III: Foster-Lyapunov criteria for continuous- 
time processes. Advances in Applied Probability , 25(3):518-548, 1993. 

[28] A. Pazy. Semigroups of Linear Operators and Applications to Partial Differential Equations . Springer- verlag 
edition, 1995. 

[29] S.T. Rachev. Probability Metrics and the Stability of Stochastic Models . John Wiley & Sons, 1991. 

[30] S. Roelly. A criterion of convergence of measure-valued processes: Application to measure branching pro- 
cesses. Stochastics , 17:43-65, 1986. 

25 



[31] S. Roelly and A. Rouault. Construction et proprietes de martingales des branchements spatiaux interactifs. 
International Statistical Review , 58(2):173-189, 1990. 

[32] G.R. Shorack and J.A. Wellner. Empirical Processes with Applications to Statistics . Wiley, New- York, 1986. 

[33] J.R. Speakman. Body size, energy metabolism and lifespan. The Journal of Experimental Biology , 208:1717- 
1730, 2005. 

[34] V.C. Tran. Large population limit and time behaviour of a stochastic particle model describing an age- 
structured population. ESAIM: PfcS , 12:345-386, 2008. 

[35] F.J.S. Wang. A central limit theorem for age- and density-dependent population processes. Stochastic 
Processes and their Applications . 5:173 193. 1977. 

[36] G.F. Webb. Theory oi Xuiiliucar iVgc-Dcpcndcnt I'opiilatioii Dynamics , volume 89 of Monographs and 
Textbooks in Pure and Applied mathematics . Marcel Dekker, inc., New York - Basel, 1985. 



26 



